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ABSTRACT 

Context. Observing the inner ejecta of a supernova is possible only in a handful of nearby supernova remnants. The core-collapse 
explosion mechanism has been extensively explored in recent models and predict large asymmetries. SN 1987A is the first modem 
stellar explosion that has been continuously observed from its beginning to the supernova remnant phase. Twenty years after the 
explosion, we are now able to observe the three-dimensional spatially resolved inner ejecta of this supernova. 

Aims. Detailed mapping of newly synthesised material and its radioactive decay daughter products sheds light on the explosion 
mechanism. This may reveal the geometry of the explosion and its connection to the equatorial ring and the outer rings around 
SN 1987A. 

Methods. We have used integral field spectroscopy to image the supernova ejecta and the equatorial ring in the emission lines of 
[Si 1] -I- [Fe 11] (A 1.64 /jm) and He 1 (A 2.058 /jm). The spectral information can be mapped into a radial velocity image revealing the 
expansion of the ejecta both as projected onto the sky and perpendicular to the sky plane. 

Results. The inner ejecta are spatially resolved in a North-South direction and are clearly asymmetric. Like the ring emission, the 
northern parts of the ejecta are blueshifted, while the material projected to the South of the supernova centre is moving away from 
us. We argue that the bulk of the ejecta is situated in the same plane as defined by the equatorial ring and does not form a bipolar 
structure as has been suggested. The exact shape of the ejecta is modelled and we find that an elongated triaxial ellipsoid fits the 
observations best. The velocity measured in the [Si 1] + [Fe II] line corresponds to ~ 3000 km s"' and is the same as the width of the 
IR [Fe II] line profiles during the first years. From our spectral analyses of the ejecta spectrum we find that most of the He I, [Si I] 
and [Fe I-II] emission originates in the core material which has undergone explosive nucleosynthesis. The He I emission may be the 
result of a-rich freeze-out if the positron energy is deposited locally. 

Conclusions. Our observations clearly indicate a non-symmetric explosion mechanism for SN 1987A. The elongation and velocity 
asymmetries point towards a large-scale spatial non-spherical distribution as predicted in recent explosion models. The orientation of 
the ejecta in the plane of the equatorial ring argues against a jet-induced explosion through the poles due to stellar rotation. 

Key words, supernovae: individual: SN 1987A - ejecta - explosions 
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1 . Introduction elusions with regard to e.g., the energy of the explosion and the 

. range of progenitors. 

— I Numerical simulations of the core collapse and explosion of , • ■ , 

massive stars have sh own that, except for progenitors with With this background it is therefore highly interesting to seek 

mass < 1OM0 diCi^ae t al. 20M), a one dimensional spheri- "^"'^h observational information about the geometry, kine- 

cally svm metric collap se do es not produce a suc cessful explo- f abundance structure of the ejecta as possible There 



X 



sion (e.g., lBuras etal.fe006l iBurrows et al.ll2007t lOessart et al.l i" '^^t^ ^as A several indications of large scale mix 



l2006h . During the last 5-10 years it has, however, become clear ^ instabilities both from t he morpho ogy and the elemert- 

that there are several effects giving rise to different kinds of mul- ta^^undance distributioii (e.g.,| Reed et al. || l995t [ Wheeler et alj 

tidimensional instabilities and convective motions. The first of ;2008|). It is howe ver, hkely that Cas A was a Type lib SN 

these is the convective motion behind the stalling shock, induced dKrause et al. || 2008[) , with only a small amount of hydrogen left, 

by the neutrino heating from the proto-neutron star (Bethe & which influences the instabilities and the general dynamic struc- 

Wilson 1985). More recendy it has been realised that a large '?'L°V.'!f-,?''T- T/'^''^^^ '^""Tl 

scale instability, usua lly known as the standing a ccretion shock °f ?N 1987 A, showing the spatial distribution of the inner ejecta 

instabiUty (SASI; .Blondin & Mezzacappal IM . is a generic P°^f radioactive decays. This is the first time that it 

feature of die core coUapse proce ss, wifli most of the power '\ P^^^^^l^ ^P^^^^^y "^f"^^ the evolution of the innermost 

in the / = 1 mode (iBlondin et alJl200lh . Several groups have ^J^^^^^ ^"^^'^^^ explosion. 

obtained successful explosions by including these multidimen- SN 1987 A is unique in being a very recent explosion, where 

sional eff'ects, although it is too eariy to draw any definite con- we have been able to follow the evolution of the SN from the 

explosion for more than 20 years. The fa ct that we know the 

* Based on observations collected at the European Southern mass of the progenitor to be around 18Mg (IWooslev et al1 l2002h 

Observatory, Chile (ESO Programme 076.D-0558) is particularly important. We have also accurate estimates of 
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the different isotopic m asses, as well as masses o f the most 
abundant elements (e.g. iFransson" & Kozmal I200I . From the 
light curve, emergence of X-rays and line profiles there were 
many indicatio ns of mixing in the ejecta (e.g. lArnett et alJ[l989t 
lMcCravl[T99l . Direct evidence for large scale instabilities came 
from HST observations, which have resolved the ejecta, and 
bo th the morp hology and the kinematics have been discussed 
by IWang et al. (2002) where it was claimed that the kinemat- 
ics and morphology as seen in the optical indicated a bipo- 
lar structure, possibly consistent with a jet outflow and possi- 
bly connected to the asym metry observed of the outer ejecta 
in early s pectro-polarimetry jCropper et al.lll988t iJeffervlfl 99 U 
IWang & Wheeler 2008i). The asymmetry in the explosion was 
suggested to have the same preferred overal l axis as detected in 
early speckle interfer ometric observations ( iMeikle et al.lll98'7t 
iNisenson et al.l [T987h and to wards the first emission knot ap- 
pearin g in the equatorial ring dSonneborn et al.iri998l: iPun et all 
I2OO2I) . The axisymmetry would include the circu mstellar ring 
structure, the outer ejecta as w ell as the inner ejecta (IWang et al.l 
I2OOI IWang & Wheeled I2OO8I) and was interpreted as possibly 
connected to a bipolar or jet-Uke structure in the explosion. In 
this paper we report ground based adaptive optics observations 
in the near-IR, which are more naturally interpreted as a prolate 
structure in the equatorial plane. 

In section|2]we describe the observations and calibrations. In 
section |3] Figures H] [5] and |7] visualise the observed spatial dis- 
tribution of the ejecta velocity field. We discuss our findings and 
compare them to previous observations of SN 1987 A in section 
m and summarise our results in section|5] 

2. Observations and Calibration 

2.1. Observations 

Integral Field Spectroscopy observations in the J, H and K band 
of the ejecta emission from SN 1987 A were obtained in Octo ber 
and November 2005 with SINFONI (lEisenhauer et al.ll2003l) on 
the Very Large Telescope (VLT) in Chile. The epochs corre- 
spond to days 6816, 6824, 6825, 6839, and 6843 since explo- 
sion. The observations were supported by Adaptive Optics (AO) 
using Star 3 as the reference source. The field of view (FOV) 
was 3" X 3" with a spatial resolution of 100x50 mas/pixel for all 
bands. 

The observations consist of 4800s in J, 4200s in H, and 5400s 
in K integrated from single exposures of 600s. All individual ex- 
posures were offset by sub-spaxel spacings to improve the spa- 
tial resolution and the object exposures were separated by sky 
exposures of equal integration time in an object-sky, sky-object 
sequence. The observations were carried out on 5 separate nights 
indicated above with airmasses between 1.4 and 1.5 for all ob- 
servations. 



2.2. Calibrations 

The data have been flat fielded, sky subtracted and com- 
bined into dat a cubes using the SINFONI pipeline version 1.2 
(ISchreiber et a l. 2004; Modigliani et al. 2007). Exposures in a 
filter for a single night (between 2 and 5 exposures) were com- 
bined using the sigma clipping routine available in the pipeline, 
where the noise is reduced by constructing the mean for each 
data point and omitting data points that were 2 sigma away from 
the mean value. 

We flux calibrated each night separately using standard stars 
(HD 76233 (B6V), CCDMJ03187 (G2V), HD 46976 (B9V), 



Table 1. The Encircled Energy (EE, here the radius) for the H & 
Kband 



EE 


50% 


80% 


H-band 


0.063" X 0.075" 


0.150" x 0.175" 


K-band 


0.075" X 0.050" 


0.150" X 0.125" 


H-band 


2.5 X 3 pixels 


6x7 pixels 


K-band 


3x2 pixels 


6x5 pixels 



HD 52447 (GOV), HD 58112 (B4V), and HD 94108 (B4V)) 
and their 2MASS values for magnitudes, types and colours 
dCohen et al.l2003l) . In the flux-calibration the observed standard 
star spectrum is divided by a template spectrum for that type of 
star in order to isolate the instrument response and telluric fea- 
tures. For the B type standard stars we superimposed hydrogen 
absorption features (scaled to the features observed in the stan- 
dard stars) on the Planck curve and used that as the stellar tem- 
plate. The Planck curves are calculated from the theoretical tem- 
peratures corresponding to the stellar type. For the G type stars 
we used the solar spectrum scaled with the observed absorption 
features in the standard star and brought it to the spectral resolu- 
tion of the standard star. 

Data from different nights were combined in order to in- 
crease the signal and to obtain a higher spatial resolution util- 
ising the sub pixel dithering. The resulting spatial resolution is 
evident in the images, where the pixel scale is 25 mas/pixel. 

We checked the pipeline wavelength calibr ation against the 
known position of the atmospheric OH lines (!Ro usselot et al.l 
|200Q) and from that determined an accuracy of the wavelength 
caUbration as the root mean square of the correction value. The 
error of the wavelength calibration is found to be 1 .25 x lO^^jum, 
2. 1 X lO^^jum and 2.4 x lO^^yum, in J, H and K, respectively. This 
corresponds to velocity uncertainties of 3 km s"', 4 km s"', and 
3 km s"'. 

The spatial resolution of our data is estabhshed using the 
point spread function (PSF) of star 3 which mostly lies within 
our field of view. The detailed knowledge of the PSF provides 
us with a measure of the residuals of the adaptive correction 
of the atmospheric disturbance and of the spatial resolution of 
the spectrograph. The PSF has an enhanced core and broad faint 
wings and we use the encircled energy metric (EE) to quantify 
the quality of the data. Since Star 3 is not completely sampled in 
the wings in all bands, we extrapolated the missing part by using 
the shape of the PSF by averaging azimuthally. We measured the 
EE in the x and y direction independently to spot any differences 
caused by the FOV image being deconstructed in the y-direction. 
The Encircled Energy (EE) for 80% and 50% of the emission 
from a point source are summarised in Table [T] The EEs are 
shown in some of the images as ellipses (akin to a beam-size) in 
the upper left corners to display the spread of the emission. 

We dereddened the spectra using the galactic extincti on law 
assum ing Ry - 3.1, and Eb-v-0.16 (Fitzpatrick & WalbornI 
Il990h for the colour excess towards SN 1987 A, based on 
Efi-v^O.lO from the LMC a nd £^-^=0.06 from the Milky Way 
dStavelev-Smith et"aL 2003h. For the re cession velocity we use 
286.7 km s ' dGroningsson et al.ll2008ah . 



3. Results 

3.1. The Ejecta Spectrum 

Our 3D AO spectroscopy allows us to acquire a spectrum of the 
ejecta and minimise the contamination from the brighter circum- 
stellar ring. Figure[T]shows the integrated spectrum of the ejecta. 
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Fig. 1. Integrated near-IR spectrum of the ejecta of SN 1987A inside a rectangular box with sides 0.575"and 0.825"encompassing 
the whole of the ejecta. The integrated area is shown in Fig. [7] 



The integrated area is a rectangular box with sides 0.575"and 
0.825"encompassing the whole of the ejecta (rectangle shown in 
Fig-EJ- In addition to the ejecta lines, there are also several nar- 
row lines coming from the equatorial ring (ER). The contamina- 
tion of the ejecta spectrum by light from the ER arises from the 
incomplete correction of the atmosphere by the adaptive optics 
which, as mentioned above, results in a PSF with broad wings. 
These narrow lines do not arise in the ejecta. 

From the ejecta we note a number of broad emission lines 
at 1.5 35yum, 1.601 yum, 1.644 //m, and 2.060//m, already identi- 
fied in lKiaer et al ] (l2007h . The 2.060/zm line is identified as He I 
(-^rest = 2.058/im). The strong lines in the H-band are blends 
of forbidden lines of singly and doubly ionised iron and singly 
ionised silicon. We also observe a weak Bry at 2.166//m. The 
region below 1.35//m is more noisy but contains several clear 
lines at 1 .20 fim, 1 .25 fj.m and 1 .27 jum. Likely identifications of 
these are discussed below. 

For the later analysis it is important to determine reliable line 
identifications, especially for the line at 1.644jum, which likely 
is a blend of [Fe II] and [Si I]. Since many of the lines originate 
from the same upper levels in Si I and Fe II, we made a simple 
model by varying the relative populations of those upper levels. 
We used Gaussian line profiles with widths of 3000 km s 
These models showed that the 1.644yum line cannot be pure 
[Fe II], since the 1.257jum and 1.322yum lines would then be 
overproduced. Rather, a mix of both [Si I] and [Fe II] is required 
for a reasonable reproduction of the spectrum. 



To obtain more specific line identifications, we have used a 
self-consistent model for th e spectral formation . This is an up- 
dated version of the code in lKozma & FranssonI (Il998ah . which 
will be discussed in detail in a forthcoming paper (Jerkstrand, 
Fransson & Kozma, 2010 in preparation). The main improve- 
ment is the addition of a Monte Carlo radiative transfer cal- 
culation, allowing for a detailed determination of the internal 
radiation field. Scattering and fluorescence are taken fully into 
account, using line opacities from NLTE solutions to the neu- 
tral and singly ionised stages of H, He, C, N, O, Ne, Na, Mg, 
Al, Si, S, Ar, Ca, Sc, Ti, V, Cr, Mn, Fe, Co and Ni. Much of 
the atomic data base is updated, both for the newly added el- 
ements (Sc, Ti, V, Cr, Mn) as well as for the previously in- 
cluded ones. Non-thermal excitations and ionisations are crucial 
at this epoch and are inclu ded for all elements along the lines of 
iKozma & FranssonI (1 19921) . 

Abundances of the di ff'erent nuclear burning zon es are taken 
from the 20 Mq model in lWooslev & Weave j (1 19951) . and we as- 
sume that these are distributed with individual filling factors in 
the core. This prescription mimics the macroscopic mixing for 
which there is abundant evidence in SN 1987A (e.g. Mc Cravl 
Il993h . At this late phase only '^"^Ti is important as an energy 
source, the positrons from which dominate the energy deposi- 
tion. We assume that these are trapped in the Fe-rich and Si-rich 
zones with 90% energy deposition in the former and 10% in the 
latter, roughly corresponding to the distribution of ^'^Ti in the ex- 
plosion model. We discuss the sensitivity to the spectrum of this 
assumption below. The total '*"'Ti mass used is 1 .0x10""* M©, and 
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Fig. 2. The upper panel shows the observed spectrum together with the total flux from our model calculation. The lower panel 
shows individual contributions together with line identifications. Dark blue is He I, red is Si I, yellow is Ca I, cyan is Fe I, green is 
Fe II. The narrow lines seen in e.g.. Pa/?, [Fe II] 1.644jum, He I 2.058 //m and Bry come from scattered light from the ring emission. 



we assume a half-life of 58.9 years (lAhmad et al.l2006l) . The ve- 
locity of the core region is 2300 km s"', in rough agreement 
with the line profiles. We do not make any assumptions about 
the composition or location of the dust, but treat it as a grey ab- 
sorber within the core. We set an optical depth of r = 1 from 
the centre to th e edge of the core, based on analyses from ear- 
lier epochs (e.g. lLucv et alJ[T991ah . We emphasise that we have 
here not attempted a full investigation of the sensitivity of the 
model to the "^^Ti mass, which depends on the assumed optical 
depth of the dust. The '^'^Ti mass should therefore only be taken 
as indicative. For a determination and extended discussion of the 
'^'^Ti mass we refer to Jerkstrand, Fransson & Kozma (2010, in 
preparation). 

In Fig. |2]we show the resulting near-IR section of the spec- 
trum. There are several interesting points to note here. First, the 
general agreement is good, given that the model has not been 
tweaked to match the data and the free parameters are few. The 
assumed '^"^Ti mass should therefore be reasonable, although, as 
discussed above, this needs to be confirmed with a full mod- 
elling of the spectrum. We also confirm that the 1 .644 yum line 
is indeed a mix of [Si I] and [Fe II], with [Si I] being the dom- 
inant component (~80%). The 1.601 jum line is essentially pure 
[Si I], 1.535 pcm. is [Fe I] and [Fe II] in similar amounts, and the 
1.50//m line is a blend of [Si I] and [Fe I]. 



The simulation also shows several [Fe I] lines. The strongest 
of these is the 1 .443 /im line, which unfortunately lies between 
the windows in which the atmosphere allows us to observe. We 
should, however, point out that the strength of the Fe I lines 
depends on uncertain contributions of non-thermal excitations. 
Collisional cross sections for these are largely lacking and are 
approximated by the Bethe formula. 

In the region below 1.35 //m we identify 1.20 yccsx with a 
blend of [Si I] 1.198,1.199,1.203 yum and [Fe I] 1.197 ^im. The 
line at 1.25 yum is mainly Fe II, and 1.27 yum is Pa/J. In general, 
this part of the spectrum is too faint in the model. The reason 
might be that the background level in the J-band is more uncer- 
tain compared to the longer wavelength bands due to the lower 
Strehl ratio. This means that the wings of the PSF, and therefore 
also the background from scattered light from the ring emission, 
is more important here than in the H- and K-bands. This could 
increase the background considerably in this band. However, the 
model also contains uncertainties (e.g. the cross sections men- 
tioned above and the positron deposition) that may also explain 
the discrepancy. 

An interesting result is that the relative strengths between 
the He I 2.058 //m line and the [Fe II] lines are well repro- 
duced. At early epochs helium from CNO nucleosynthesis gives 
rise to the 2.058 //m emission. In our model most of the he- 
lium and iron emission comes from the iron core, although we 
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Fig.3. Left panel: image of the 1.644|im [Si I]+[Fe II] feature. The white dashed ellipse indicates the apparent shape of the inner 
ring, which is centred on (0,0) marked with a cross. The ejecta shape is indicated by the black dashed ellipse. The ellipses in the 
top left corner show the 50% (80% in red) encircled energy area from a point source. The colour bar gives the flux intensity in 
10"'^ erg s"' cm"^. The right panel shows the line profile of the 1.644yL(m feature extracted at three different positions shown in the 
left panel. The blue curve corresponds to the upper most extraction box, the green to the middle box, and the red curve to the bottom 
box. 



cannot exclude some contribution by helium from the hydro- 
gen burning region, mixed into the core. Evidence for this was 
seen at earlier epochs from the line profile of the He I lines 
jKozma & FranssonI 1998b|), as w ell as in 2D hydrodynamic 
simulations ( Kifonidis et al.l 120061) . The He in the iron core is 
produced by the g-rich fre eze-out. According to the models by 
IWooslev & Weaveii(ll995l) the He abundance in the iron core can 
be ~5Q% by number One may in this context note that the He 
I 2.058 fim emission is substantially stronger than the Bry line 
(Fig.[Tli, which is mainly driven by the recombination freeze-out, 
unless positrons leak into this region. The optical depth to the 
gamma-rays is small, and only a small fraction of the gamma- 
ray energy is absorbed by the ejecta. Recombination freeze-out 
effects also affect the helium rich gas from the envelope, and add 
to the He I emission from the core. In our model most of the He I 
emission is, however, produced by the Fe core. A more quantita- 
tive assessment of the He I emission requires a detailed spectral 
calculation. 

It is not clear whether the positrons are absorbed by the zones 
containing the "^^Ti, where they were emitted, or i f they propa- 
gate t o other zones before being absorbed (see e.g. lChugai et al.l 
Il997h . In order to take this uncertainty into account, we explored 
the consequences of this assumption by computing a model, 
where we assumed that the positrons were deposited in propor- 
tion to the total number of electrons in each zone, rather than at 
the same location as they were emitted. 

This did not result in any major changes in the model spec- 
trum, although the origin of the emission in some case were dif- 
ferent. The 1.601 i-im and 1.644 fim lines were still dominated 
by Si I in this model, however now produced by the silicon in 
the oxygen zones. This alternative model yielded a lower flux 
for these lines by a factor ~ 2. The He I 2.058 /jm line flux did 
not change much between the two models, but the source of the 
emission changed from the Fe/He zone to the He and H zones in 
the alternative model. 



3.2. The Ejecta Geometry 

We focus in the following on the 1.64-4-fim and He I 2.058 //m 
lines in our spatial analysis of the velocity field as, from our 
data, these lines provide the strongest constraints on the nucle- 
osynthesis. 

First we investigate the apparent shape of the ejecta. The 
left panel of Fig. [3] shows an image of the brightest ejecta line 
at 1.644yum. The small cross is positioned at the co-ordinates 
(0,0), which corresponds to cf = 05h 35m 28.105s and 6 - - 
69° 16' 10.99"(J2000.0, with the astrometry calibrat ed with re- 
spect to HST images and the coordinates given by IWest et all 
ll987l:IWalbornetalJl993h . This is the centre of the white dashed 
line, that follows the shape of the inner ring. The shape of the 
ejecta is elliptical and we find from the isophotes a ratio of the 
major to minor axis of 1 .8 + 0. 17 (see black dashed ellipse in left 
panel of Fig. [3]). The axis of symmetry is at a position angle of 
15° + 0.9°, which agre es quite well \ yith th e 14° + 5° found from 
HST observations by IWang et aH (l2002h and ISugerman et al.l 
(I2005h . 

With 3D spectroscopy we have spectral information for each 
pixel in the image in the left panel of Fig. [3] The right panel 
of Fig. |3] shows the spectra at three different positions indicated 
by the squares in the left panel of Fig. [3] The extraction is per- 
formed using an averaging box function 0.1 arcseconds on the 
side which maximises the signal to noise while not significantly 
degrading the resolution. The different spatial positions of the 
spectra lead to different line profiles of the 1.644yum line. The 
southernmost pixel (red curve) has most red-shifted emission, 
the northernmost pixel (blue curve) has most blue-shifted emis- 
sion, and the middle pixel (green curve) has both red and blue- 
shifted emission and thus a very broad profile. The spectra of 
the North and South pixels alone indicate a bipolar structure of 
the elongated ejecta. The middle spectrum supports this with its 
width. 
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Fig. 4. Images of the 1.644yL(m line of SN 1987 A and its circumstellar ring in different velocity bins along the line of sight, which 
show the spatial distribution of the different velocities for that line. The fluxes are in units of 10""* erg s ' cirT^. The colour bar to 
the left indicates the scaling for the middle bin and the colourbar to the right the remaining bins. The contours trace the intensity 
levels, spaced by 5 x 10 erg s"' cm"^. The asterisk marks the centre of the inner ring (as in Fig. O. 



In order to use the spectral information to derive ejecta ge- 
ometry along the line of sight we briefly summarise the expected 
kinematic structure of freely expanding ejecta. Because of the 
kinematic nature of the explosion, once pressure effects have 
become unimportant, the expansion is expected to be homol- 
ogous. With a homologous velocity field, i.e., V{r) oc r, con- 
stant velocities along the line of sight are surfaces perpendicular 
to this direction. The images at the different wavelengths of an 
emission line therefore represent a tomography of the structure 
of the ejecta along the line of sight. In particular, for a spheri- 
cally symmetric ejecta we expect the images at a given veloc- 
ity along the line of sight, V~, to be concentric disks with radius 
p - Roi^-iVz/Vo)^y^^, where Vo is the velocity at the maximum 
radius Rq. Departures from spherical symmetry will show up as 
shifts in intensity in the centre of the images and a non-circular 
form. 

To investigate the asymmetry of the ejecta we show in Figure 
ID images of the 1.644jum line of SN 1987A in different veloc- 
ity intervals. Each frame in the figure shows a separate part of 
the ejecta line profile. The divisions are made in velocity, where 
the line is integrated over 1000 km s The upper panels show 
the blue-shifted emission and the lower panels show the red- 
shifted emission. The contours show the intensity levels spaced 
by 5 X 10"'^ erg s"' cm"^. 

The integrated image of He I 2.058 jum clearly shows elon- 
gated ejecta similar to that in the 1.644jum line. We have com- 



pared the 1.644jL/m and He I images by taking the ratio of the 
two, but do not find any significant difference in extent or inten- 
sity distribution. 

Although not the subject of this paper (see instead lKiaer et alJ 

l200l . we note the dominance of the ring in the lowest 
velocity bins of both lines. Furthermore, we also note the 
stronger emission from the southern part in the red-shifted 
[500, 1500] km s"' bin, while the northern part dominates the 
blue-shifted [-1500,500] km s"' bin. Although the lines from 
the shocked ring are of i ntermediate width (200-400 km s"') 
dGroningsson et alll2008ah . we see that the ring is visible in all 
the velocity bins. This is caused by the continuum and other lines 
from the ring and/or a broadening of them due to the shock in- 
teraction in the ring. The most important result, which is new, is, 
however, that in all velocity bins there is a clear shift in the posi- 
tion of the maximum intensity between the blue and red sides of 
the ejecta, consistent with the velocity pattern seen for the emis- 
sion from the ring. This illustrates in a nice way the power of 
AO-supported integral field spectroscopy. 

We can also study the geometry in other lines. The spatial 
distribution of the kinematics for the He 1 2.058 fim line is shown 
in Figure|5] Although the emission for the ejecta is much fainter 
and the ring is brighter in He I 2.058 yum than in the 1.644jum 
line, we observe the same velocity pattern in the emission aris- 
ing from the ejecta, with the blue-shifted emission appearing to 
the North and the red-shifted component to the South. The ring 
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Fig. 5. Images of the He I 2.058 /im Une of SN 1987A and its circumstellar ring in different velocity bins along the line of sight, 
which show the spatial distribution of the different velocities for that line. The units are the same as in Fig.H] The contours show 
the intensity levels, spaced by 2.5 x 10"^ erg s ' cm"^, and white are saturated areas (to enhance the ejecta against the bright ring). 



emission seen in the [-1500,-2 500] km s ' bin is most likely 
due to the [Fe II] 2.046/im line (iKiaer et al.ll2"007l) . 

3.3. Simulations of the ejecta kinematics 

In order to understand the 3D shape of the ejecta using the spa- 
tial and velocity information in Fig. |4] we constructed a model 
based on a homogeneous ellipsoid. This is obviously a simplifi- 
cation in view of e.g., the possible patchy du st obscuration, see n 
as a 'hole' in the HST images of the ejecta jWang et alj|2002h . 
and also compared to the comp l ex structure seen in re alistic sim- 
ulations dKifonidis et al.ll2006t [Hammer et al.ll2010l) . As a first 
approximation to a non-spherical distribution the homogeneous 
ellipsoid provides the necessary constraints for a qualitative de- 
scription of the ejecta geometry. We have allowed the axis ratios 
and the Eulerian angles of the ellipsoid to vary. Finally, the im- 
ages have been convolved with the observed PSF profile for Star 
3 in the H-band. 

Fig.|6]shows the spatial distribution of the emission in differ- 
ent velocity bins, scaled to represent that of the observations in 
Fig. ID The centroid of the emission is determined by a combina- 
tion of the shape of the ejecta, as represented by the axis ratios, 
and the orientation in space of the emission region. Observations 
during the first years showed strong evidence for dust in the 
ejecta (iLucv et al.l [l991bt IWooden et al.l Il993h . Dust will ob- 
scure emission from the far side of the eject (the red part of the 
integrated along the line of sight emission). However, since we 
spatially and spectrally resolve the emission, dust could only im- 



pact quantitative measures of the fluxes of lines arising from the 
ejecta. The main conclusion of our study, namely, that the ejecta 
are confined primarily to the equatorial plane, is unlikely to be 
altered by the efi'ects of dust obscuration. We therefore do not 
include this effect in the simulations. 

The best fit to the observations is achieved with a triaxial 
body with axis ratios x : y : z = 3 : 2 : 1, with Eulerian angles 
6 = 90° (inclination), (p = 25° (line of nodes) and i// = 80° (angle 
of major axis). Here w e follow the convention for the angles by 
iGoldstein etaP (l2002h . This corresponds to a position angle of 
10° and a tilt out of the plane of the sky by 25°. The position 
angle is consistent within the errors with that determined directly 
from the image. 

Comparing the model (Fig.|6]l with the observations we see 
that there is a general qualitative agreement between the two as 
seen in Fig.|4l with the centroid of emission indeed shifting from 
north to south between the blue and red. Also the axis ratio be- 
tween the projected images on the plane of the sky agrees well. 
Quantitatively there are, however, differences as expected for 
this simplified model. While the red-shifted velocity bins com- 
pare reasonably well, the blue-shifted bins have a smaller spatial 
extent, corresponding to an ellipsoid, where the major axis has 
a larger inclination from the plane of the sky (i.e. viewed more 
'head on'). We discuss the ejecta geometry more in section|4] 

This comparative analysis suggests that the ejecta shape does 
not follow a simple ellipsoid. Rather, that the shape of the blue- 
shifted part is different from the shape of the red-shifted part, or 
at least that the viewing angles of the two are not the same. 
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Fig. 7. Images of the spectral and spatial distribution of the ejecta Unes. Top left: The [Si I] / [Fe II] 1.644jum line, integrated over 
the spectral range shown in the spectrum to the right. The encircled energy for 50% and 80% of the emission from a star is indicated 
by the ellipses. See section|2]for details. The colour bar is in units of 10"'^ erg s"' cm"^, and the contours are [20, 40, 60, 80, 100, 
120, 140] in the same unit. The white box indicates the integration area for the spectrum shown to the right. Right: The spectrum, 
with the bottom panel keeping the spatial position along the North-South axis, with respect to the asterisk in the left image. A 
projection of the encircled energy is shown for 50% and 80%. The colour scale follows the colour bar to the right with the unit 
10^'^ erg s~' cm"^ A"' □" ', and the contours are [2, 5, 8, 11, 14, .., 41] in the same unit. The integrated spectrum for the whole 
area is displayed above, where the flux is in the unit 10"'^ erg s"' cm"^ A"' . Bottom: Same as above but for the He I 2.058 jum 
line. 




Fig. 6. Simulations of the ejecta intensity from a uniformly 
emitting ellipsoid for the same velocity intervals as for Fig. |4] 
Negative velocities are in the top row and positive in the bottom. 
The cross marks the centre of the ejecta. See text for details of 
the model. 



3.4. The Ejecta Kinematics 

The elongated ejecta visible in images of SN 1987A show that 
the North part of the ejecta is predominately blue-shifted and 
the South part is predominately red-shifted. Figure [7] shows the 
spectral and spatial distribution of the [Si I]H-[Fe II] 1.644//m 
line and the He I 2.058 fim line. The left panels of the figure 



are image maps integrated over the spectral range from the right 
panels, where the top image is the 1.644jum line and the bottom 
image the He I 2.058 jum line. The right panels show the line 
profiles for the given lines and underneath the spatial distribu- 
tion of the emission line profile along the North-South axis, also 
referred to as the spectral image. The images are centred on the 
asterisk at coordinates (0,0), which marks the centre of ellipti- 
cal appearance of the inner ring. The spectral images follow this 
convention for the positional zero point. The Encircled Energy 
(EE) for 80% and 50% of the emission from a point source are 
indicated by the ellipses (see Section|2]for details). 

The projection of the EE ellipses are displayed in the spectral 
image, keeping in mind that rather than displaying a positional 
error, they display the distribution of the emission. The white 
boxes indicate the integration area for the spectrum shown in 
the right panels. The area of the boxes are chosen in order to 
minimise the contamination from the ring emission while opti- 
mising the emission from the ejecta. The spectral image retains 
the North-South axis, showing us the spatial distribution of the 
red and blue shifted material. Any information about the spatial 
distribution along the East-West axis is lost in this integration. 
Reversing the direction of extraction such that the information 
in the EW direction is preserved did not result in any additional 
kinematic information. 

The similarity of the spectral image, and thus the kinemat- 
ics, for the 1.644yum and He I lines suggests that He is inter- 
mixed with the Fe/Si. This is consistent with an origin of the 
helium from the a-rich freeze-out in the iron core, as we dis- 
cussed in connection to the spectral analysis (section 13.11 ). As 
was described there, we can, however, not exclude that mixing 
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Fig. 8. Spectral images of the fainter ejecta lines in the H-band, 
showing the spatial distribution along the North-South axis. Top; 
the [Fe II]/[Fe I] 1.53 yum line. Bottom: the [Fe II]+[Si I] 1.60 yum 
line. The colours follow the colour bars to the right in units of 
10"'^ erg s ' cm"^ A ' □" ' and the contours trace the intensity 
levels, spaced by 3 xlO"'** erg s"' cm"^ A"' The spatial 

distribution clearly follows the one of the bright 1.644/^m line 
in Fig. It] 



of the helium from hydrogen burning may give the same spatial 
distribution as the iron core. 

Fig. [8] shows the spectral images of two fainter ejecta lines 
in the H-band at 1.53 //m and 1.60//m. Here we see that these 
lines show the same spatial distribution of the velocities along 
the North-South line (y-axis) as the 1.644/im and He I lines. 

Note that we only observe gas which is heated and excited by 
either radioactivity or the emission from the shock. The latter is 
likely to be important mainly for the outer high velocity ejecta. 
The inner ejecta we are discussing in this paper is most likely 
mainly excited by the radioactive material. At these late epochs 
the ejecta is transparent to the gamma-rays, and the positrons 
from the "^'*Ti decay dominate the excitation. Because these are 
likely to deposit most of their energy locally the emission seen 
from the inner ejecta should mainly reflect the distribution of 
radioactive "^^Ti. The explosion imprinted a velocity structure on 
the whole ejecta, and it is this velocity structure that we now see 
illuminated from within by radioactive decay. 



4. Discussion 

The observed kinematics and geometry of the ejecta of 
SN 1987 A at late times clearly indicate that the emission has 
an orientation generally similar to the equatorial ring. Fig. |9] 
shows the position of the ejecta emission along the line of sight 
and compared to the position of the equatorial ring. We find the 
inner ejecta of SN 1987A to be distinctly aspherical, with the 
blue-shifted material displaced to the North and the red-shifted 
material predominately luminous to the South (Figs.|4]|5]|7]|8]l. 
The spatial extent of the red-shifted (southern) material is signif- 
icantly larger than that of the blue-shifted material. If we, in spite 
of the differences between the blue and red shifted material, ap- 
proximate the shape of the whole ejecta with that of an ellipsoid, 
we find that the orientation of the major axis is PA=15° and with 
a tilt out of the plane of the sky of ~ 25°. 

The maximum velocities observed are about 2500 km s"' to 
nearly 3000 km s '(Fig. |7]l, corresponding quite well with the 
velocities inferred from the line shapes of the radioacti ve decay 
lines and the observed Fe lines in early spectra (e.g. iMcCravl 
We observe material, which is still heated by radioactiv- 
ity, hence must come from a region mixed with elements syn- 



Fig. 9. Schematic view of the ejecta distribution of the ejecta rel- 
ative to the the ring as seen in the 1.644//m line. Top: The dis- 
tances from the centre are given in meters, the observer is located 
on the left. This figure demonstrates that the ejecta mostly lie in 
the same plane as defined by the equatorial ring. The squares 
only give approximate emission centres. In reality the emission 
is more diffuse, see text for details. 



thesised in the explosion. With the inner ejecta not perpendicu- 
lar to the equatorial ring, we find a con figuration, which d evi- 
ates from the bipolar model proposed bv lWang et al.l (l2002h and 
Wang & Wheeler (20081) . 



Wang et al] ( l2002h discussed the geometry and kinematics of 



the ejecta based on HST observations obtained in August 1999, 
some six years bef ore our SINFONI observations. The HST 
image discussed by IWang et al] (l2002h was taken through the 
F439W filter and is dominated by emission in [Fe I] and [Fe II]. 
It should therefore be directly comparable to our [Si I]-i-[Fe II] 
SINFONI image. The position angle of the ejecta was found to 
be 14° + 5°, consistent with our results. Wang et al. (2002) also 
discuss the kinematics based on the line profile of the [Ca II] 
A 7300 line. By applying 0.1 "slits along the North-South axis, 
they found that the peak velocity of the northern part was close 
to zero, while the southern ejecta has a positive velocity of 
~ 1700 km s"'. These numbers should be compared to our av- 
erage velocities of -1400 km s ' and ~700 km s~', respectively 
(cf. Figl3]i. However, as noted bv lWang et al.l (l2002h . there may 
be systematic errors in their velocities, either from the uncer- 
tainty in the rest wavelength of the line ([Ca II] A 7300 vs. [O II] 
A 7320), in the wavelength calibration, or the exact positioning 
of the slit. Considering these uncertainties, our measurements 
may th erefore be consistent, // the wavelengths by IWang et aP 
(l2002h are blueshifted by -1000 km s ' . We note, however, that 
our spectra have higher S/N, and that we can trace the line pro- 
files more reliably. In addition, the ejecta has expanded in the 
seven years between the two observations (the supernova age in- 
creased by a factor of nearly 1.5) and hence a better separation 
of the different emission regions can be achieved. Finally, thanks 
to the integral field spectroscopy, we can image the ejecta emis- 
sion over the whole ejecta, as well as in specific lines, which 
represent a considerable advantage compared to the HST filter 
observations. 

The observed geometry of the inner ejecta strongly sup- 
ports a large-scale instability, like SASI, in the explosion. The 
strong mixing of material and the asymmetry of the explosion in 
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SN 1987 A, as already indicated by the bolometric light curve, 
the Bochum event, line profiles and the early emergence of j- 
and X-rays is here clearly confirmed. The kinematics of the in- 
ner ejecta show that the explosion has defined an orientation, 
which is close to the equatorial plane of the progenitor star as de- 
fined by the circumstellar rings. This would invalidate jets pro- 
duced through the poles, if the equatorial ring really defines a 
rotation axis of the progenitor star. 

In addition, the inner ejecta appear not to be located in a 
plane, but rather we see two extensions at different radial ve- 
locities and projected spatial extent. The strongest emission in 
the northern extension is located at about -1400 km s^'and at 
0. 15"from the explosion centre, while the southern part displays 
the strongest emission at a velocity of ~500 km s^'but extends 
out to nearly 0.3". Fig. |9] displays the situation of the ejecta 
relative to the equatorial ring. The ejecta velocities have been 
converted to the line-of-sight positions at the time of the ob- 
servations (6840 days) and the angular distance in declination 
((5DEC) assuming a distance to SN 1987A of 50 kpc. The ring 
position was determined from mea sured offset from the centre 
and assuming an inclination of 41° ( ISugerman e t al. 2005). This 
figure should be seen as a more detailed 2D complement to Fig. 
|6] which, however, gives a full 3-D representation, although with 
the approximation of an ellipsoid. The general orientation of the 
major axis is consistent within the errors of these two represen- 
tations. 

Since the inner ejecta seem not to be co-planar, we may also 
speculate about the direction of the kick received by the neu- 
tron star. Presumably, the inner ejecta are the part, which impart 
the momentum of the large hydrodynamic instabilities onto the 
proto-neutron star The connection between the SASl instabili- 
ties on scales of a fewxlO^ km and the large s cale asymmetry 
we see is n ot obvious. The first calculations bv iKifonidis et alj 
(l2006l) and lHammer etal.l (1201 Ol) show that the SASI is impor- 
tant for triggering of the instabilities, although it is not clear how 
much of the initial instability survives the dynamics of the explo- 
sion on the much longer time scale of hours and days until the 
hydrodynamic structure is frozen-in. It is therefore too early to 
speculate about a specific direction for the kick of the neutron 
star. It should be recalled that the images reflect the distribution 
of the positron input and although the positrons are stipulated to 
have a limited travelling range their distribution cannot be taken 
to be the distribution of mass in the core. 

As was discussed in Sect. 13.31 dust affects the observed flux 
from especially the red side of the emission. The limited signal- 
to-noise and the fact that most lines are blended make a detailed 
analysis of the dust difficult. The fact that we do observe sub- 
stantial emission from red-shifted velocities (e.g., Figs.[3]and|7ll 
shows that the dust can only block the core partially at most, as 
it has earlier in the evolution. Even so, the line profiles of the 
rather clean [Si I] 1.60 and 1.64 jum lines show clear deficien- 
cies of the red sides. We have not attempted any investigation 
to see whether these asymmetries have quantitatively changed 
compared to earlier epochs (e.g. iFassia et al . 2002). mainly be- 
cause of the limited resolution and signal-to-noise of these ob- 
servations. 

5. Conclusions 

The outer ejecta of SN 1987A travelling at around 30000 km s"' 
have reached the outer ring a decade ago, and have become 
visible mostly in the r everse shocks generated by the collision 
with the ring material (McCravl 20031 Groningsson et"ani2006i 
l2008blal : lKisr et al.ll2007HZanardo et al.ll2010h . The inner ejecta 



are still heated b y the radioactive decays (m ostly positrons from 
the "*"*Ti decavs: iFransson & Kozmall2002h and have expanded 
enough to be spatially resolved. The extent of this material cor- 
responds to about 3000 km s"' . Our spectral analysis shows that 
most of the helium emission comes from the helium produced 
in the a-rich freeze-out in the core. This is also supported by the 
similar spatial distribution of the He I and [Si I] + [Fe II] images. 
Helium zone material mixed into the core may, however, also 
contribute to the He I emission. 

We have confirmed the asymmetric shape of the inner ejecta 
in SN 1987 A and have shown that it is confined roughly to a sim- 
ilar plane as the equatorial ring. Also, the northern and southern 
lobes are not symmetric and show slightly different radial ve- 
locities, which points towards two separate angles in the line of 
sight of the emission sites. Both thes e arguments are ag ainst a 
jet-induced explosion as favoured bv lWang et al.l (l2002h due to 
an explosion propagating through the poles of a rotating star. 
Instead, the shape of the inner ejecta are fully consistent with 
what is expected from the large instabilities predicted in recent 
explosion models of core-collapse supernovae. The SINFONI 
observations of SN 1987 A can be seen as a direct observational 
confirmation of these models. 
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